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Abstract. The two-channel Anderson lattice model which has SU{2) ® 5(7(2) symmetry is of relevance 
to understanding of the magnetic, quadrupolar and superconducting phases in Ui_2,Tha;Bei3 or Pr base 
skutterudite compounds such as PrFe4Pi2 or PrOs4Sbi2. Possible unconventional superconducting phases 
of the model are explored. They are characterized by a composite order parameter comprising of a local 
magnetic or quadrupolar moment and a triplet conduction electron Cooper-pair. This binding of local 
degrees of freedom removes the entropy of the non Fermi-liquid normal state. We find superconducting 
transitions in the intermediate valence regime which are suppressed in the stable moment regime. The gap 
function is non analytic and odd in frequency: a pseudo-gap develops in the conduction electron density 
of states which vanishes as close to u — 0. In the strong intermediate valent regime, the gap function 
acquires an additional k-dependence. 

PACS. 74.20.-z Theories and models of superconducting state - 74.20.Mn Non-conventional mechanisms 
- 75.10.Lp Band and itinerant models - 75. 30. Mb Valence fluctuation, Kondo lattice, and heavy-fermion 
phenomena 



1 Introduction 

Superconductivity or superfluidity of ■^He [0 is understood 
to be a thermodynamic phase in which spin 1/2 Fermions 
are subject to Cooper-pair correlations and condense. The 
Fermions are either electrons or ^He atoms. Superconduc- 
tivity is well described by the Bardeen-Cooper-Schrieffer 
(BCS) theory |^ when the fluctuations of the order param- 
eter are small. The Eliashberg theory of superconductivity 
is applicable in superconductors when high and low 
energy scales are well separated. Contrary to the BCS the- 
ory, this theory accounts for the retardation of the interac- 
tion mediating Bosons. On a large time scale, the electron 
repulsion is effectively overcome in standard superconduc- 
tors by the electron-phonon interaction. The small param- 
eter needed for the applicability of the Eliashberg theory 
is given by the ratio of the Debye frequency and the Fermi 
energy. It is well known that Eliashberg theory cannot be 
used in the '^He problem p| since no separation of energy 
scales can be found here. He atoms pair in spin triplets 
and orbital momentum of L = 1 in contrast to standard 
superconductors which usually develop spin singlets pair- 
ing. The rich phase diagram of ^He is directly related to 
this anisotropic pair-condensation. Pairs with L = 1 have 
nodes in space and, therefore, can effectively overcome the 
hard-core repulsion and gain condensation energy. 

Heavy Fermion materials have drawn much atten- 
tion since the discovery of superconductivity in CeCu2Si2 



Q , which seems to be characterized by an anisotropic or- 
der parameter with a symmetry yet to be determined. 
Another example for an anisotropic order parameter are 
the high temperature cuprate superconductors. Their cur- 
rent carrying subsystem consists of layered copper ox- 
ide perovskite structures. Strong experimental evidence 
has been compiled for an unconventional order param- 
eter in these materials, with /2 symmetry of the form 
cos{kj;a) — cos{kya) or (fc^ — ky) which is usually called 
"d-wave" superconductivity by many authors. A super- 
conductor is called conventional if the order parameter 
has the full symmetry of the Fermi-surface; the order pa- 
rameter still might be anisotropic but it transforms in 
accordance with the trivial irreducible representation of 
the point-group A. Otherwise, the order parameter is 
called unconventional. Recently, the superconductivity of 
Sr2Ru04 has generated strong interest due to structural 
similarities between this compound and the High-Tc ma- 
terials 0. It has been suggested that Sr2Ru04 might be 
an example of an unconventional spin triplet superconduc- 
tor 1^,^ in contrast to the High-Tc materials where a /2 
symmetry is favoured. In spite of the progress in our un- 
derstanding of Heavy Fermion materials, there has been 
little success in determining the symmetry of their order 
parameters in the superconducting phase. Huth et. al. re- 
cently reported experimental evidence which indicates an 
electronically mediated pair-condensation mechanism in 
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UPd2Al3 Our theoretical understanding of unconven- 
tional superconductivity is still far from complete in spite 
of more than 20 years of research. Theories were developed 
along two basic lines: either the superconducting phase of 
microscopic models was investigated using, at that time, 
state of the art approximations in order to reveal possi- 
ble pairing mechanisms, or phenomenological approaches, 
based on the Ginsburg-Landau theory, were constructed 
in order to give new insight into the possible symmetries of 
the order parameter and their impact on the measurable 
quantities pl| . 

Two channel Anderson and Kondo lattice models ex- 
hibit non- Fermi liquid behaviour in the paramagnetic phase 
drive n by unquenched and fluctuating local degrees of free- 
dom This phase is characterized by a large resid- 
ual resistivity and entropy, and ill defined electronic quasi- 
particles. Fermi liquid physics is restored by cooperative 
ordering or apphed magnetic field or stress In par- 
ticular, for a non Kramer's doublet crystal field state in 
U^^ or Pr'^"'" ions, a two channel Kondo effect is possible 
]l5t . The magnetic "spin" of the electrons serves as chan- 
nel index in this case. UBcis and PrFePi2 are prominent 
candidates for a quadrupolar Kondo lattice description, as 
their enormous resistivity (> 100/iJ7cm) is removed only 
by phase transitions (superconductivity in UBcia, anti- 
ferroquadrupolar order in PrFePi2). Indeed, it has been 
shown that commensurate and incommensurate orbital 
ordering, as well as ferromagnetism are possible in this 
model depending upon coupling strength and filling 
Evidence of a first order transition to an odd frequency 
pairing state in the Kondo limit (near integral valence) 
has been adduced, which is a singlet in both spin and 
channel indices, and no q value was preferred for the cen- 
ter of mass momentum 



qvi 



Our paper is organized as follows. In section 0, we re- 
view the two-channel Anderson lattice model This 
model is rotational invariant in spin and channel space 
and, therefore, has SU{2) ^ SU{2) symmetry. We use es- 
tablished standard methods developed by others ||l^,f3|, 

p,ra,|T7ira,pi,po|,i2iii2~ 



to obtain an approximate solu- 
tion of its one- and two-particle properties. In section ||, 
all possible spin and channel symmetry sectors of Cooper- 
pairs are discussed. Two of the four sectors require odd- 
frequency gap functions if we restrict ourselves to It- 
dependence with 7^1 symmetry. We discuss the concept of 
a composite order parameter in section 5.3. A connection 



between the composite order para met er and the Cooper- 
pair susceptibility is made in Sec. 3.4. The statements of 
Section ^ are exact and independent of the any approx- 
imation. The solution of the homogeneous phase of the 
model is used to obtain the superconducting transi- 
tion temperature defined by the divergence of the pair- 
susceptibility which will be the subject of Sec. ^. In sec- 
tion H we present the self-consistency equations for the 
triplet/triplet 8x8 Nambu matrix Green function, which 
is the so-called dynamical mean field theory (DMFT) p^ , 
analog of the Eliashberg equation. However, no 
Migdal theorem is required since the full local irreducible 
pairing vertex is used. This non-perturbative set of equa- 



tions lead to quasi-particle and anomalous Green func- 
tions which are gapped at cj = 0, as shown in Sec. 5.4. The 
gap function is non-analytic in frequency, and, therefore, 
Heid's theorem is not applicable. This superconduct- 
ing phase belongs to a minimum of the free energy. 



2 Theory 

2.1 Model 

Uranium based Heavy Fermion materials exhibit a lot of 
interesting and unusual physical properties: coexistence 
of magnetism with superconductivity and non-Fermi liq- 
uid behaviour in the paramagnetic phase. The Uranium 
5f-shell states in Heavy Fermion compounds can be 5/^ 
(U^+), 5/2 (U'*+) and 5/^ {\J^+). Experimentally, how- 
ever, 

U5+ 

is energetically unlikely. The Hund's rules ground 
state multiplets are given by J = 5/2, J — A and J — 9/2 
respectively and split in a cubic symmetry ||2^ accord- 
ing to the irreducible representationss of the point group. 
If we consider only valence fiuctuations between the low- 
est irreducible representationss of 5/^ ^ 5/^ or 5/^ <-^- 
5/^, the Uranium 5/-shell is modelled by a magnetic dou- 
blet r(i{r^) above a non-Kramers doublet ground state 
which has a quadrupolar, and hence orbital character. 
Only conduction electrons with Ig symmetry, labelled by 
a spin a and a locally defined orbital index a, couple to 
these ionic states according to a simple group-theoretical 
selection rule • The same arguments also apply 

to Pr^+ (J = 4) and Pr^^ (J = 5/2) in a cubic crys- 
talline environment which is found in the skutterudite ma- 
terials PrFe4Pi2, which shows anitferro-quadrupolar order 
around 2K, and PrOs4Sbi2 where a superconducting tran- 
sition was found at 2K |2^ . If we restrict ourselves to lo- 
cal hybridization and coupling to two generate band - an 
assumption which we will justify within the DMFT the 
lattice model on a cubic lattice 

kcTQ (71/ OLV 

+ ;^E(^"4a.^r3-.,A.+/..c) (1) 

^ VQi(7 

has an SU{2) (g) SU{2) symmetric with respect the spin cr 
and channel a space. ^, are the usual Hubbard projec- 
tion operators ||2^ at lattice site f which describe transi- 
tions from local 4(5)/ atomic state 7' to the state 7. 

The Hamiltonian (|l|) contains three different energy 
scales: the band- width D of the conduction bands £aa, 
the inter-configurational energy Saa = E„ — Ea and the 
so-called Anderson or hybridization width A = irpoV^. po 
is the conduction electron density of states at the chemical 
potential = 0. If not otherwise stated, we will use A as 
the energy-scale throughout this paper. 

2.2 Dynamical Mean Field Theory 

The fundamental assumption of the DMFT is the locality 
of single particle self-energy which allows the mapping of a 
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complicated lattice problem onto an effective site problem 
|p^,po| which is embedded in a fictitious bath of conduc- 
tion electrons. This bath, the so-called media, takes into 
account the angular averaged but energy dependent influ- 
ence of the lattice onto this effective site. The DMFT self- 
consistency condition (SCC) demands the equality of the 
k-summed lattice Green function with the Green function 
of the effective site. Then, only one-particle irreducible 
loops, which describe the hopping of a single electron from 
and to the effective site, enters the SCC. Any two parti- 
cle irreducible loops - correlated hopping of two electrons 
from and back onto the site - are factorized in a prod- 
uct of one-particle loops: higher order correlations which 
also contribute to a local self-energy, are neglected. This 
assumption becomes exact in the limit of infinity spatial 
dimensions |l],|2|] . 

Even if we had taken into account a band dispersion 
£kaaa' , uou-diagonal in the orbital index a, the mapping 
onto an effective site requires two locally defined, degen- 
erate effective conduction bands which are symmetry pro- 
jected linear combinations of the Wannier state of the con- 
duction bands. Cox has shown that only those hybridize 
with the local /-states on symmetry grounds Since 
the local dynamics is governed by the correlations induced 
by the /-states, it is consistent with the DMFT to restrict 
ourselves to diagonal band dispersions Eka-a in Eqn. (|I|). 
We use the so-called Non-Crossing Approximation (NCA) 
|^,|2[|3^ as a solver for the effective impurity through- 
out the paper since it has been established as standard 
technique for solving the two-channel lattice model 

E- 

2.2.1 Single-Particle Green Function 

The conduction electron Green function Gaa (k, z) for the 
model obeys the exact equation of motion: 



which connects it to the /-Green function 

F„.(k,2) = ^E^"''''' «X'U.\K.-c.» (^) 



^ ^kacT 

(2) 



Since FQ,o-(k, z) can always be written as 



IT4 



•2 S^kacT 



(3) 



(4) 



where in DMFT F-^^{k, z) = F-^{z) = z-ec.a- S[^{z), 
the DMFT solution for the single-particle Green function 
of two-channel periodic Anderson model is determined by 
the solution of 



1 



(5) 



-E 
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Fig. 1. Bethe-Salpeter equation for the two-particle Green 
function IJ^^ , F denotes the irreducible vertex. 



The local self-energy IJl^^{z) is calculated by a SIAM with 
an effective conduction-electron density of states p{oj) = 
'^mA{Lo~i5)/\TiV^\. Therefore, the self-energy of the con- 
duction electrons 17^^ (k, z) = z — ekatr — G~^(k, z) is given 



by 



1 -|- GTaa-{z) Z — Eaa ^ ^aa{z) 



(6) 



as a consequence of the equation of motion (H) where 
Taa{z) — \ya^Faa{z). Equs. (||-||) are valid for any num- 
ber of channels. 

In contrast to the single channel Kondo effect, the 
ground state cannot be mapped onto a renormalized Fermi 
liquid in the two channel case. Here, the effective coupling 
constant flows to an intermediate value. The correspond- 
ing fixed point has non-Fermi liquid properties and is un- 
stable against a magnetic and a channel field . 

The single particle and transport properties of this 
model [p^ show non-Fermi liquid-behaviour within the 
DMFT approximation which can be seen immediately by 
the analytic form of the conduction electron self-energy 
(^: independent on the method which is used to solve 
the effective impurity, the local /-self-energy and, 
therefore, also self-energy of the conduction electrons, must 
have a finite imaginary part even for T — > in the para- 
magnetic phase since the local T-matrix is pinned to half 
of the unitary limit [^|3^,^. The scattering originates 
in the two-channel Kondo effect in the impurity model 

: the local quadrupolar moments of the /s states are 
dynamically screened by the conduction electrons orbital 
moments, if < E^. In the case of E^ > E^ orbital 
and spin degrees of freedom are interchanged, and the lo- 
cal magnetic moments are screened instead. The predicted 
excitation gap at low temperatures in the optical conduc- 
tivity has been found in UBcis ^6| . 

This result of a finite resistivity in the paramagnetic 
phase, originating in magnetic scattering at the /-sites due 
to the absents of singlet ground state, indicates the insta- 
bility of this phase towards those phase transitions which 
will remove the residual entropy of the lattice. We showed 
in a previous publication [l4|| , that the model has a rich 
phase diagram of spin-density, quadrupolar-density, and 
anti-ferromagnetic and anti-quadrupolar order for larger 
|eQo-| > A- For the intermediate valent regime {\£aa\ < ^) 
ferromagnetism is found for small band fillings. We sus- 
pect, that superconductivity dominates the intermediate 
valent regime. 
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2.2.2 Particle-Particle Green Function 



where 



The main focus of this investigation is devoted to the 
study of superconducting transitions. One indicator of a 
phase transition is a divergence of the static order param- 
eter susceptibility. Therefore, we are interested in Fourier 
components of space-time susceptibilities of the type 



Xo.,oAt) = - <T{OS)Ol)> , 



(7) 



where Oi will be a symmetry adapted Cooper pair op- 
erator CiaaCia'a-'- The lattice susceptibilities are obtained 
as 



Xo(q, i^) 



1 



x7T(icj„,k; ja;,„,k')(q, I/) , 



FoiiLUn,k)F^{iu;,n,k')i8) 



where Fo{iWn,k) is the Fourier component of the struc- 
ture factor of the operator Oi{t) and lUiujn, k; iujrm k')(q, v) 
denotes the two-particle Green function. The summation 
over spin and channel indices has been omitted in the 
following secti on f or clarity. It will be introduced and dis- 
cussed in Sec. 3.4. 



The two-particle Green function 77 obeys a Bethe- 
Salpeter equation, diagrammatically depicted in Fig. |l]. 
The lines are fully renormalized one-particle Green func- 
tions. A scaling analysis has proven that the total irre- 
ducible vertex has to bepurely local in the limit of infinite 
spatial dimensions |^,^. Total irreducibility means that 
each diagrammatic part remains connected when any two 
electron lines are cut. Although, we use the textbook def- 
inition of two-particle irreducibility with respect to Feyn- 
man perturbation theory, the structure of the equations 
are invariant of the underlying perturbation theory since 
the Bethe-Salpeter equation is the two-particle equivalent 
to the Dyson equation. 

It was noticed early in the development of the DMFT 
p^,^ that only the local part of the irreducible two- 
particle vertex enters the Bethe-Salpeter equation for the 
two-particle Green function. Moreover, the renormaliza- 
tion of the one-particle self-energy by the two-particle 
vertex is already taken into account in DMFT ||2l|,^. 
This statement, however, is only true inside the homoge- 
neous (paramagnetic) phase. The susceptibility becomes 
divergent near the phase transition temperature Tc where 
at least one irreducible vertex diverges and becomes k- 
dependent. The scaling argument fails in this case. There- 
fore, the self-energy E acquires an additional, possible k- 
dependent, contribution in the symmetry broken phase. 
This has been widely used to investigate possible mag- 
netic and superconduction phases in various models |^, 

With a local irreducible pp-vertex 7^ in the homoge- 
neous phase, the particle-particle Green function 7J reads 



7J(ia;„,k;itj,„,k')(q,w„) = 

n,m 



(9) 



-|-7T/ree(k, q, iuJn-, iVn)I{iLUn, i^m, ^n, Q) 

x7T/ree(k',q, , (10) 



77/ree(k,q, itJ„, W„) = Gk+q(it^„ +iVn)G-]^{iuJn) (11) 

and the reducible vertex I(iujn, i^m, ii^m Q ) obeys the Bethe- 
Salpeter equation 

I{iU!n, iuJrn, i^n, Q) = r{iuJn, il^,n\ i^n) (12) 

xxo(jw„/; w„,q)7(iw„/,iw„)(w„,q) . 
Herein, the particle-particle bubble xoii^n] ii^n, Q ) is given 

by 

Xo(«w„;ii/„,q) = -i- ^ 7T/ree(k, q, zcj„, w„) . (13) 

k 

Energy w„ and momentum qof the center of mass (COM) 
are conserved. It was shown p3 that the COM-momcntum 
only enters through a scalar ?7(q) 



1 



(14) 



on a hyper-cubic lattice in d dimensions. In this case, 
Eqn (02) is casted into a matrix equation in Matsubara- 
frequency space 

£(ii^„,q) =^(«i^„) + -^£(jM„)x^(w„,q)X(w„,q) , (15) 



which is formally solved by 



I(*'^n,q) 



1 " ^£(*'^«kg(«'^"'q) 



Liiiyn) ■ (16) 



The DMFT mapping of the lattice problem onto an 
effective impurity problem is used in order to calculate 
the irreducible vertex r^^{ivn)- In formal equivalence to 
Eqn. (p^), the equation 



- +^,„^^£(-«)S,„,(-n) , (17) 

depicted in Fig. |l|, holds for the local two-particle Green 
function of the effective site. It has been shown that the 
irreducible vertices -r(w„) in ( p^ and in ( p6| ) are identical. 
They may be obtained by matrix inversion in Matsubara 
frequency space 



p — =loc '■o 



(18) 



once the free local pp-Green function 

Xioc(iw„, iuJm\ Wn) = G(iw„ -I- iVn)G{-iuJn)kur„i^^ (19) 

and the local pp-Green function 77, (wn) are known. 

^=loc ^ ^ 
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aa — iLO„ 



a a ium + iv-n 




Fig. 2. Local particle-particle Green function for the two- 
channel model in NCA plotted on an imaginary time axis, 
from bottom to top. The wiggly vertical line represents the 
Fz states, the double straight the magnetic /e states and the 
horizonal lines the hybridizating conduction electrons. 



2.2.3 Local Two-Particle Green Function 

The local two-particle /-Green function Uf^iuJm ioJm] ii^n) 
for the periodic Anderson model has been calculated within 
the NCA in the limit of two local valence configurations, 
one with N electrons and one with A'' -1-1 electrons [Q. For 
the two-channel model, the diagrammatic contribution is 
shown in Fig. |^ and reads: 

Zeff Jc 

xP^a' {z + iuJn - ii^m)Pcr' (z - iuJm) (20) 

where the contour circumvents all singularities of the ionic 
propagators P-y{z). Zeff is the partition function of the 
effective site. The indices on the bottom of 71/ label the 
spin and channel of the incoming while the indices on the 
top label the outgoing particle-hole pair. Since there is 
one external line crossings when 77j?^ is inserted into the 
Bethe-Salpeter Eqn. ([l7|), the overall sign must be nega- 
tive. 

In general, Eqn. ( pO| ) has to be numerically evaluated. 
At r = 0, however, this integral can be solved analyticly 
for \iujn\ <C A for the two-channel Anderson model if we 
use the ionic propagators assuming that the lattice 
self-consistency does not change the low-frequency thresh- 
old behaviour. This is justified in the paramagnetic phase 
for an energy independent DOS close to w = 0. A rather 
complicated calculation yields 

rr f ■ n\ ^ ( ^ V sign(t^n) • sign(w™) 



niax(w„, 
min(Lj„,cj„) 



where F{^,x < 1) is the elliptic integral of the first kind. 
For low frequencies iujn and icom, 77/(«w„, zti;,„; 0) turns 
out to be independent of the low energy scale Tk- The con- 
duction electrons over-screen the local spin: a compensa- 
tion of the spin can never be achieved, and the conduction 
electrons continue to be incoherently scattered for T ^ 0. 
The two-particle scattering is governed by low lying exci- 
tations independent of the cross-over energy scale Tk, and 
the temperature T appears to be the only relevant energy 
scale. This approximation is valid for |i^|,T < Tk- 

Since the large Coulomb repulsion prevents the ionic 
shell from accepting or donating two electrons at the same 
time, a local pair correlation at site i can only be gener- 
ated sequentially. The first electron has to leave the /- 
shell before another conduction electron can enter it, see 
Fig. 1^. Therefore, the particle-particle /-Green function 
tends to have nodes as a function of time. 7T/(iCi;„, iwn', 0) 
becomes an odd function of both of its frequency argu- 
ments iwn and icom in the two-channel model due to the 
two exchangeable doublets. In other models with locally 
restricted charge fluctuations, we expect that odd and 
even components with respect to frequency will be present. 

The periodic Anderson model contains two types of 
electrons: localized and itinerate ones. Since the bands are 
non-interacting in the absence of hybridization, particle- 
particle correlations only occurs through hybridization with 
localized states in ionic shells. We can express the conduction- 
band particle-particle Green function in DMFT as 

loc ~r Xc,loc TlocXc, loc ! 

(22) 

where energy arguments have been dropped on the r.h.s. 
The local two-particle T-matrix Ti^^^ cumulant is defined 

by 

Tioc{iuJn,it^,n\iVn) = (U f {iuJri, iuJn, il^n) 

-l3F{-iujn)F{iiyn + it^„)5„,m)(23) 

Xc,loc 

the media particle-particle Green function is 

Xc,loc = G{-iuJn)G{iuJn + iVn) ■ 

rif factorizes in a product of two /-Green functions in 
the absence of local interactions, and the two-particle T- 
matrix cumulant vanishes as expected. Using Eqn. (|l^ ) 
and (^ ) we obtain 



-r =A-T, A— 
13='' =/3='°==l- 

A^x , [x , ] 



1 



(24) 



(25) 



The matrix A measures the degree of renormalization of 
the band electrons. It is of order 0(1) at high temperatures 
and becomes moderately enhanced in the low temperature 
regime. Hence, the irreducible vertex F ~ T, .If we 

^ ^ =c ^loc 

define a renormalized T-matrix 21;^^ = A H^^^ A, the 
irreducible vertex assumes the same structure as the one- 



particle self-energy I]aa{z) (^ 



max(w„, LOri 



(21) 



F =T, 



(26) 
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Sector 


Spin 


Channel 


Frequency 


# Orderp. 


I (SsCs) 


singlet 


singlet 


odd 


1 


II (SsCt) 


singlet 


triplet 


even 


3 


III (StCc) 


triplet 


singlet 


even 


3 


IV (StCt) 


triplet 


triplet 


odd 


9 



Table 1. Classification of the 16 possible order parameters 
for local Cooper pairs. Pauli's principle determines the fre- 
quency behaviour of the symmetrized anomalous Green func- 
tion F{z). Odd frequency implies {F{z) = ~F{~z), even fre- 
quency F{z) = F{ — z). Only Cooper pairs with ^(k) trans- 
forming according to the irreducible representations Fi are 
taken into account. Sector I and IV show odd frequency be- 
haviour, while in sector II and III conventional even frequency 
pairing is found. 



where q denotes the momentum of the Cooper pair. 

Since 7 contains the spin and orbital degrees aa, the 
possible symmetry properties of the pair susceptibility can 
be assigned to the four sectors listed in table ^ These are 
associated with certain subspaces of the product space of 
spin and channel degrees of freedoms. Each of them are 
spanned by a product of singlet and triplet states as in- 
dicated. Since the total symmetry under exchange of the 
electrons has to be odd due to Pauli's principle, the fre- 
quency dependence of the symmetrized anomalous Green 
function F{q, iuJn) 



F{q, iujn) 



77'k 



(31) 



3 Superconductivity 

3.1 Classification of the Cooper-Pair Operators 

It is assumed that only the electrons in the conduction 
electron sub-system are directly involved in the formation 
of Cooper pairs in the condensed phase. The quantum 
numbers 7 and 7' label their spin and orbital degrees of 
freedom. If a Cooper pair has an internal structure mod- 
elled by a structure factor S'(2r), its pair operator reads 



P^,^' (R) = ^ S'(2r)c^ (R - y)cy (R + r) 



(27) 



has to be odd in sector I and IV; hereby, A^^' denotes the 
symmetry form-factor of the appropriate sector. 

The most general pair operator in k-space in sector I 
(spin singlet/channel singlet) is given by 

P^^(k, k') = S{^^) aaCi,o.aC-W-a-a , (32) 



from which a symmetrized pair operator of the form (|2 
can be obtained by 



P^^(q)=^P^'^(k + q/2,k-q/2). 



(33) 



r = (ri — r2)/2 determines the location of one electron 
relative to the center of mass R=(ri-)-r2)/2. The trans- 
formed pair operator P-y^^'(q) 

P-y.f' (q) = E (k - q/2)cy (-k - q/2) (28) 

k 

describes a pair in reciprocal space with a center of mass 
momentum q. 

Since in DMFT the irreducible particle-particle (pp) 
vertex P is also a q-independent quantity, compare Eqns. (j|) 
to (p^), the pair susceptibility 

xp,7,7',7",7"'(q,^) - (r(P7,7'(q)W4',7"'(q))) ' (29) 

is only influenced by the local interaction if the structure 
factor S'(k) transforms with Pi, the trivial representation 
of the point group. Hence, solely superconducting insta- 
bilities with a conventional order parameter which also 
transforms according to Pi can be obtained within DMFT 
through (|29|). This includes, however, Cooper pair form- 
factors such as 5(k) cx £k, sometimes called extended s- 
wave, which will play the leading role in the extreme mixed 
valent regime where Eaa = 0. Therefore, we will ignore 
unconventional spatial form factor symmetries which, of 
course, cannot be ruled out in the model. 

The corresponding one-particle anomalous Green func- 
tion P-y_^/(k, q, z) is given by the Fourier transform of 



When we introduce the transposed bi-spinor operator 
'0"^(k) = (ck+ti Ck+x, Ck-|, Ck-x); the general pair opera- 
tor in sector II (spin singlet/channel triplet) can be writ- 
ten as 

P^^*(k,k') 5(^^)^^(k)z^^ *Z^ZV^(-k') (34) 



and as 

P*^^(k, k') = S{^±^W(y)^Z^ «^,^^(-k') (35) 

in sector III (spin triplet/channel singlet), where g_ acts 
in the spin sector and t in the channel sector. The nine 
pair operators of sector IV ( spin and channel triplet) are 
given by the matrix with the elements 

i^*-*(k,k') = S{^±^)V{^^)^^^^^ *X,X,^(~k') . (36) 



3.2 Introduction to Odd-Frequency Pairing 



^'7,7'(k, q,T) 



-(7'(c7,k-q/2(-r)c^ 



-k-q/2 



)> 



(30) 



In section 2.2.3, it was shown that the irreducible vertex P 
is proportional to the local /-particle-particle propagator 
Flf in leading order, and that /7j is odd in the incoming 
and outgoing frequencies at low temperatures, Eqn (pl|). 
Hence, we expect to find a tendency towards so-called odd- 
frequency superconductivity in model (|l|). 

Berezinskii has pointed out in the context of ■^He that 
s-wave/spin triplet or p- wave/spin singlet pairing, both 
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having a temporal node, are also compatible with Pauli's 
principle it only requires a total parity of -1. The 
anomalous Green function has to be an odd function, 
with respect to frequency, in the above mentioned cases. 
However, it quickly became apparent that ^He orders in 
p-wave (L=l), spin triplet (S=l) pairs which show even 
frequency behaviour. Berezinskii's idea was forgotten, but 
was later revitalized by Balatsky et. al. in the context 
of superconducting phases of the Hubbard-model. Similar 
to magnetic phases, a superconduction ordered phase can 
have an order parameter whose modulation in real space 
is described by a finite wave vector q. This wave vector 
corresponds to a finite center of mass momentum q of the 
Cooper pairs. Coleman et. al. investigated the pos- 
sibility of an odd-frequency solution for a Kondo lattice 
model and found a stable solution for such a staggered or- 
der parameter within a mean-field theory. This solution, 
however, breaks spin-rotational invariance. The order pa- 
rameter consists of a local spin contribution bound to a 
conduction electron and transforms as an 5 = 1/2 object. 
Abrahams et. al. argued that the order parameter of 
an odd-frequency superconductor has to be a composite 
object, i. e. a Cooper pair bound to a localized spin, in or- 
der to produce a positive Meissner effect. Therefore, time 
reversal symmetry is not broken by odd-frequency pairing 
since both components of the order parameter, the Cooper 
pair and the spin operator, have negative parity with re- 
spect to the time-reversal operator. Held, however, showed 
that odd-frequency superconductors are thermodynami- 
cally instable against q- modulation . He assumed that 
the normal-phase was described by Fermi-liquid theory 
and that the order parameter did not have singularities. 
Emery and Kivelson |35j found an enhancement of a com- 
posite pair operator susceptibility using the solution of 
the two-channel Kondo-model at the Toulouse point. This 
could indicate a tendency towards an odd-frequency phase 
in a corresponding lattice model. 

Since the mechanism of superconductivity in strongly 
correlated electron systems is still subject to an ongo- 
ing debate, it is of great importance to explore this un- 
usual concept in particular since, as will be seen below, it 
is favoured by orbital degrees of freedom. Whether odd- 
frequency pairing is actually realized in nature, is not yet 
clear. 



3.3 Composite Order Parameter and its Physical 
Interpretation 

Pauli's principle demands that superconductivity in the 
spin singlet, channel singlet (SsCs) sector or in the spin 
triplet, channel triplet (StCt) with a Ji form-factor must 
have an odd anomalous Green function as indicated in 
table ^ Assuming pairing with q = 0, the local pair ex- 
pectation value 

{Pss{(l = 0)) = ^J2Y1 aaSik){c^aaC-i,-c.-a) = 
k aa 

(37) 



vanishes with the anti-commutator {ckaa, c-k-a-u} — 0, 
which reflects Pauli's principle. Hence, this expectation 
value cannot serve as the superconducting order param- 
eter of section I and IV. This is a peculiarity of odd- 
frequency superconductivity. 

Abrahams et. al. noticed that the time derivative D^^i — 
{-^Cj{t)cj') is non- vanishing in the odd frequency super- 
conducting phase and may possibly be furnished as an 
order parameter Since the time derivative operator 
does not change the symmetry of the spin and orbital part, 
the order parameter can still be classified according to the 
symmetry of the conduction band Cooper pairs listed in 
table H This new kind of order parameter is formed by 
linear combinations of -D^^', which will be called O be- 
low. The derivative -^Cj{t) is equivalent to the commu- 
tator [H,Cj{t)]. In the case of the two-channel periodic 
Anderson model, ■j^Cj{T)cy consists of a product of two 
operators: a conduction electron pair operator and a lo- 
cal magnetic or channel spin operator both generated by 
the hybridization term between the band electrons and 
the local /-states. Thus, one arrives at a composite order 
parameter, which correlates local and itinerant degrees of 
freedom: formation of band Cooper pairs is stimulated in 
the presence of a local pair resonating between definite 
spin and channel states with proper symmetry. 

Since there is only one order parameter component Oss 
in sector I, it must be a global singlet. The only possibil- 
ity to generate a singlet from a vector operator such as 
the local /-spin S is to form a scalar product involving 
another vector operator in the spin space. It has to be the 
conduction electron pair operator P*''' given by Eqn. (|3^ ) 
which is a scalar in orbital and a vector in spin space: 

Oss = sign{AE) ((S P*-^) - (r P^^*)) . (38) 

The second term in this expression, the first term's sym- 
metrical counterpart, combines a local /-shell orbital spin 
T and a pair operator P**'* which is a vector in orbital 
space and a scalar in spin space to form a global singlet. 
Both terms are expected to contribute since the Hamilto- 
nian (|l|) models local fluctuations between spin and orbital 
doublets at each lattice site. In this way, the "pair-wave 
function" in sector I realizes correlations between the dy- 
namics of the local spin and a conduction electron pair, 
which are compliant with Pauli's principle. 

The dimensionless order parameter Oss is invariant 
under particle-hole transformation and simultaneous ex- 
change of the two local doublets. This is also a fundamen- 
tal symmetry of the Hamiltonian. The exchange of spin 
and channel spin is compensated by the sign change of 
AE — Ea- — Ea = £aa being independent of spin and 
channel in the absence of a magnetic field or lattice dis- 
tortion. Hence, the superconducting order parameter is 
not dependent on whether the magnetic or the quadrupo- 
lar doublet has the lower energy. Moreover, the composite 
order parameter meets the objective that the supercon- 
ducting state must bind local magnetic or orbital degrees 
of freedom in order to quench the residual entropy present 
in the two-channel model. 
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The pair operators which enter Oss in Eqn. (jS^) trans- 
form according to the irreducible representationss of sec- 
tors II and III. This reflects the fact that Pauh's principle 
allows that their expectation value may stay finite. Of 
course, the symmetry of the composite order parameter 
is determined by the scalar product as a whole and not 
by the symmetry of the individual operators; therefore, 
the order parameter transforms according to sector I. A 
simple mean field decoupling of Oss , as given by 



O,, -sign(Z\i;)((S)(P* 



(39) 



would essentially treat the local and itinerant correlations 
as independent quantities and would have to be inter- 
preted as a mixture of phases of sectors II and III rather 
then superconductivity in sector I. Correlations between 
the conduction band pairs and the local spin and orbital 
degrees of freedom are essential for generating the unusual 
phase in sector I. Nevertheless, this decoupling might still 
be a way to understand the experimentally observed co- 
existence of magnetism and superconductivity in certain 
Uranium Heavy Fermion compounds. An antiferromag- 
netic spin order, in conjunction with superconductivity 
with a zero center of mass momentum of the Cooper pair, 
would translate into a modulated order parameter Oss, 
with wave vector q given by the magnetic modulation 
vector. Alternatively, one could speculate about a finite 
momentum Cooper-pair [^5| bound to a modulated spin 
or orbital ordered state in such a way that the total mo- 
mentum of the order-parameter is zero. 
With the definitions 



(40) 



= ]^ E e^""'*^"^ko.G,,_,,_,^,^,(*c.„) ,(41) 

it is straightforward to derive exact relations for the sym- 
metrized order parameter for the singlet/singlet channel 
in the form 



Os 



V 



using equations of motion as sketched in appendix ^ The 
order parameter has two contributions: the term Dss de- 
scribes an isotropic local Cooper pair, i. e. S'(k) = 1, while 
Tss is formed by nearest neighbour pairs, i. e. S'(k) = 
Ek/i, sometimes called extended s-wave pairing. The It- 
summation on the r.h.s. of (|l]) leaves T^^,' finite smce 
^cka„,c^ka'a' (*^") ^ '5'(k) for the case of extended s-wave 
pairing. 

The composite order parameter Oss contains no nat- 
ural distinction between isotropic and anisotropic pairs 
with a form-factor eu/t. It is interesting to note that for 
the intermediate valence regime \E^—Ea\ — > 0, the nearest 
neighbour Cooper pairs with extended s-wave symmetry 
dominate since the first term on the r.h.s. of (143) vanishes 



for \AE\ — > 0. This can be interpreted as condensation 
of delocalized pairs. There will not only be nodes in fre- 
quency in this regime but also nodes in k-space. 

An exact equation similar to (42) can be derived for 
the tensor order parameter in the triplet-triplet sector 



O 



= sign(Zii?) {{s,P;^') - lj,Pl 



\AE\ 

1/2 



1,3 



sign{AE) 
y2 



T 



(43) 
(44) 



where i,j — x,y,z, i{j) indexes the spin(channel) space, 
Pj'* or P-*'" is a component of the singlet /triplet or triplet /singlet 
pairs (HJ, ^ and the symmetrized isotropic (Dij) and 
nearest neighbour pair (Tij) expectation values are ob- 
tained by the symmetrization rules ( ^6|) in combination 
with (|4^) and (|l|): 



(45) 



3.4 Composite Order Parameter and 
Pair-Susceptibility 

As explained above, the order parameter is formed by a 
composite operator combining a conduction band pair op- 
erator in singlet/triplet channel with a local spin or chan- 
nel operator in the two odd-frequency sectors I and II. It 
measures the correlation between a local spin and a con- 
duction electron pair. To verify the existence of such an 
unconventional phase in the model, it is necessary to ei- 
ther show that there is a self-consistent solution of the 
order parameter equation in the symmetry broken phase 
or search for divergencies of the order parameter suscep- 
tibility xo(q) — Go.o^{iVn — 0, q). It is cumbersome, 
however, to obtain a tractable expression for xo{^.) since 
it involves a three-particle Green function due to the com- 
posite nature of O. On the other hand, Eqns. ( ^2| ) and (Q) 
connect the expectation value of the order parameter to 
the derivative of the symmetrized pair operator ^cjj,(T)cy 

which takes the form P°!^f{iujn) = iuinCj{iuin)Cj'{—iLLjn) in 
the Matsubara frequency space. Therefore, it is equivalent 
to investigate the pair susceptibilities 



xp(q) - ^ 



E 



G podd(i^^-) p\ add(i^^){ivn — 0,q) 



(46) 

They are obtained from two-particle Green functions whose 
four external lines are connected to a frequency dependent 
and, in the case of "extended s-wave" k-dependent vertex, 
generated by P°^f{iujn) and P\°f''^ {iujm) ■ 

Up to now, all statements in section ^ are exact and 
independent of any approximation. In the following, how- 
ever, we will calculate the two-particle Green functions 
only within the DMFT. Knowing the irreducible particle- 
particle vertex P for the appropriate symmetry, and only 
considering isotropic pairs, the pair susceptibility reads in 
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general 



Xp(q,iM„) = - 



E 



iUni-i^m) X(q,«l^n) (47) 



1 



Since r_ is local within the DMFT, only gap functions 
transforming according to Fi symmetry can be investi- 
gated. All other symmetries would require a microscopic 
solution of the model beyond the DMFT or the restric- 
tion to simple BCS mean-field theories. Especially latter 
approaches almost always confirms the ansatz and the as- 
sumptions made at the beginning of the calculations, and, 
therefore, it is not clear to us whether real physical insight 
is gained. 

The phase transition temperature Tc and the frequency 
dependence of the anomalous self-energy (l>{iLUn) is ob- 
tained by solving the eigenvalue equation 



M|A,M) = A|A,M) , 
where the symmetric matrix M is defined as 



M=^/A(q, 0)^00), /x(q,0) 



(48) 



(49) 



The divergent contribution to the pair susceptibility at 
ivn = 0, i. e. 



xp(q) = 



1 



1- A, 



-regular terms 



(50) 



P(q) 



1. Hereby, cx„ 



stems from the largest eigenvalue A„ 

given by (A|^x(q,0)|w„) and |u;„) = (• • • , w„, • • (see 

appendix |b| for more details). This is the DMFT-analog to 
the Eliashberg equation in the standard theory of super- 
conductivity and was first used by Jarrell and co-workers 



4 Results for the Phase Boundary 

As discussed before, all two-particle Green functions are 
evaluated within DMFT in the paramagnetic phase using 
standard procedures The temperature Tc at which 

the pair susceptibility diverges defines the boundary be- 
tween the paramagnetic and the superconducting phase, 
if it is the highest instability temperature. 



4.1 Sector I: Spin Singlet/Channel Singlet (SsCs) 

The pair operator in the odd-frequency SsCs sector 
2 



(q,zi.„) = — ^ azc.„5(k) (51) 




T/A 



Fig. 3. Contribution pSsCs ^j^g ^^^^ susceptibility Xf>''^"{<i) 
from the dominant maximal eigenvalue |A| vs temperature in 
the singlet/singlet (SsCs) sector for q = and q = 7r(l, 1, 1) 
(upper plot). The lower plot shows the total pair susceptibility 
Xp given by Eqn (p7|). Parameters: Saa- ~ —3 and ric — 2. 



is derived from (p^), where S{k) — 1 for isotropic and 
S'(k) = Su/V of "extended s-wave" pairs. Hence, the pair 
susceptibility xf? consists of a spin diagonal and a spin or 
channel flip contribution. Now, the spin a nd cha nnel in- 

in order 



dices, which have been neglected in section 2.2.2 



to focus on the structure rather than on details, are of im- 
portance: the resulting coupled Bethe-Salpeter equations 
yield 



x|,'"^"(q) = iUniLUrn[n'^'''{iuj„,iLUm;q,0) 



/52 

-i7'='=(iw„,iu;™;q,0)] 



(52) 



in the absence of a magnetic field and isotropic Cooper- 
pairs. IJ'^''^ describes spin and channel diagonal particle- 
particle propagation; 77*^^ the exchange of spin or channel 
index. The coupled equations are solved for the particle- 
particle Green function TT'** = jjdir _ jjex fj.Qj^ which the 
effective pp- vertex Ff' entering (^7|-^9|) 



=loc. 



j-jdir YI^^ 

^=loc ^=loc 



(53) 
(54) 



is obtained for the singlet/singlet sector, where Gc denotes 
the conduction electron Green function. 

Jarrell et. al. investigated the pair susceptibility for 
the two-channel Kondo model in the strong coupling limit 
using QMC fl^ . They did not find a maximum positive 
eigenvalue A™ with Am(T'c) = 1 for a given finite tem- 
perature for q = or q = 7r(l,l,l) which is equivalent 
to ?7(q) = — 1 (staggered pairing). However, they observed 
that a single eigenvalue A approaches — cx3. By lowering the 
temperature further, the associated eigenvector, i. e. the 
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symmetry form-factor, corresponds to a very large posi- 
tive eigenvalue which decreases with decreasing tempera- 
ture leading to a sign change of the contribution PcsSc to 
■^SsCs |-ggg pjg 3 j^g£ Since an eigenvalue A > 1 

was found, it was interpreted as an indicator for a possible 
superconducting instability using a linearized Eliashberg 
equation. The simultaneous sign change of the staggered 
and the q = contribution is clearly visible in Fig. 3 
of Ref. ||l3|] which looks very similar to the upper plot 
in Fig. M. Since A appears to jump suddenly to a value 
A > 1, it has been interpreted as a locally driven first 
order phase transition . The susceptibility correspond- 
ing to the composite order parameter Xp^'~^^ j which would 
indicate a true instability of the whole system, was not 
investigated. 

In general, the local conduction-electron particle-particle 
Green functions are obtained from a local two-particle t- 
matrix equation analogous to the ph-Green function ( p2|) 
for models with non-interacting conduction bands. The lo- 
cal connected part of the two-particle T-matrix cumulant 
is given by 



(7Tf'^(it^„,ic^™;0) 



n,m } 



.;0) 

(55) 



Within NCA, the spin or channel exchanging two-particle 
f-Green function HJ^ , depicted in Fig. ^, is of the order of 

0{v^ = 1), where v — V p. However, the local two-particle 
f-Green function, diagonal in spin and channel indices, 
turns out to be of the order O(w^) and, therefore, is ne- 
glected leading to a total contribution of 77^^(iw„, iujm \ 0) 
given by . We solved the Eliashberg- like e qua tion (|8|) 
and confirmed the results of Jarrell et. al. |l^ for the 
two-channel periodic Anderson model. The temperature 
dependence of p-5sCs^ depicted in the upper plot of Fig. 
^, shows excellent agreement with the published results 
for the two-channel Kondo model (see Fig. 3 of Ref. jl^). 

In Appendix we show that the sign change of a con- 
tribution to the pair susceptibility is equivalent to find- 
ing a non-trivial solution to the linearized DMFT self- 
consistency condition (SCC) in the s ymme try broken phase. 
This will be discussed later in Sec |3.1.l| in more detail. 
However, the form factor which was found for the Kondo 
model [ |l3[ and the two-channel periodic Anderson model 
in our calculation appears to be singular at w = 0: F(uj) cx 
1/uj. This has the effect that a linearization of the DMFT 
equations, on which the linearized Eliashberg equation is 
based on, is not justified for |w| 0, and the full SCC 
has to be considered before any conclusion about the ex- 
istence of a symmetry broken state can be drawn for sec- 
tor I. Additionally, the susceptibility of the composite or- 
der parameter xf.*'^'* is three orders of magnitude larger 
than the particular component P^'^'^'^ exhibiting the sign 
change which, therefore, might be an irrelevant contribu- 
tion to the total pair-susceptibility of the composite order 
parameter Oss- 

We investigated the SsCs sector for a large range of 
parameters and did not detect a phase transition or a 
sign change in Xp'''^*- The sign change of the contribution 

pScCs 

to Xp^*^^ has been viewed as indication of a first 




T/A 



Fig. 4. Contribution P''"-'* to the pair-susceptibility xp 
from the largest positive eigenvalue vs temperature in the 
triplet/triplet (StCt) sector for q = 7r(l, 1, 1) (staggered) and 
q = (upper plot), and total pair-susceptibility xp, Eqn. (^^, 
vs temperature (lower plot). Parameters: Eaa = —3 and ric = 2. 



order phase transition by Jarrell and coworkers . How- 
ever, we always find a higher transition temperature in the 
spin-triplet, channel triplet sector for the investigated pa- 
rameters of the two-channel Anderson lattice model. We 
believe that the charge fluctuations on the /-shell plays a 
major role in the formation of the superconducting phase 
of the Anderson lattice model as can be seen from Eqns. 
( ^2|) and ( ^ ) while these charge fiuctuations are absent in 
the Kondo-lattice model. The superconducting gap func- 
tion is dominated by an so-called "extended s-wave" ( Pi) 
form factor S'(k) = J^a cos{ka) parameter in the extreme 
intermediate valence regime, AE = 0. 



4.2 Sector II (SsCt) and Sector III (StCs) 

We looked into the even frequency eigenspace of Ff^ (|5 
and P" ~ 



-P" = E 



f3 



E 



(56) 

which belong to the StCs and the SsCt sector, respec- 
tively. E^_^^Jn,m = {5n,7n +'5m^„)/2 projccts out the odd- 
frequency components. No instability was found in either 
sector. This was expected from the analytic form of two- 
particle T-matrix ( |2l| ) which is as ympto tically odd in its 
frequency arguments (see also Sec. 2.2.3). 



4.3 Sector IV: Spin Triplet/Channel Triplet (StCt) 

The spin triplet, channel triplet sector, StCt, is character- 
ized by a tensor order parameter with nine components 
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0.004 
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1.5 



2.5 



(b) 



Fig. 5. Reduced transition temperature tc — Tc/T* (a) as a 
function of the band filling Uc for four different values of \eacr\ 
and (b) the dependence of Tc and T* on ric for e^o- = — 3Z\. 
The components of superconducting order parameter Oij are 
invariant under Sacr —Sacr and a simultaneously performed 
particle-hole transformation, as can be clearly seen in the plot. 
Even though Tc increases for decreasing |£ao-|, the ratio Tc/T* 
decreases since increasing charge fluctuation tends to destabi- 
lize the superconductivity. 



( ^ ) (see also table Q). Some of the two-particle Green 
functions in this sector involve both a direct and an ex- 
change contribution yielding a matrix representation in 
Matsubara frequency space which shows an odd-frequency 
dependence in order to obey Pauli's principle, for example: 
the spin and channel diagonal Green function. This struc- 
ture is already revealed in the free local particle-particle 
Green function 

X^o"(lW„, icj„i; iVn) = Gaa{i>^n + il'n)G „a{-i>^n) 

for ivn = 0. Consequently, the matrix inversion is mathe- 
matically ill-defined in the particle-particle version of Eqn. 
( p6| ) since the determinant of an odd matrix in Matsubara 
frequency space must vanish. This problem can be circum- 



vented by defining a matrix inversion for odd-frequency 
matrices in relation to an odd frequency unitary matrix 
Eodd{n,m) = {Sn,m - <5n,~m)/2. This analysis, however, 
is only necessary for some of the nine different pp-Green 
functions which are needed for the diagonal order pa- 
rameter susceptibilities Xp'i^* — Gq (0). By careful 

analysis of all triplet two-particle Green functions, it can 
be shown that all odd-frequency pair-susceptibilities are 
equal in DMFT, and the relevant irreducible vertex is 
given by 



E 



^=loc 



jjex 
^=loc 



E 

^=.odd 



(58) 

where E_^^^ projects out the odd-frequency components. 

For the two-channel Kondo-lattice model, Xp^^z 
investigated by Jarrell et. al. p^ , and no superconduct- 
ing instability was found for Kondo couplings of J/t > 0.4. 
A different picture emerges in the two-channel periodic 
Anderson model. In Fig. ^, the total pair-susceptibility 
y^stct (flower curves) and the contribution stemming from 
the largest eigenvalue P'S'*C'* (upper curves) are shown in 
the stable moment regime for half-filling and for the cen- 
ter of mass momenta q = and q = 7r(l, 1, 1). A sign 
change is observed in both quantities at a temperature of 
T « 8 • 10~^Z\ and for q = 0, which originates from the 
divergency of Pstct as Amax 1- Consequently, we find a 
clear tendency towards a second order phase transition in 
the triplet/triplet sector for an uniform order parameter. 
The pair-susceptibility of a staggered pair is much smaller 
than x|.*'^*(q — 0) and stays continuous in the vicinity of 
Tc- The qualitative difference of the Kondo-lattice model 
and periodic Anderson model is the possibility of change 
fluctuation on the /-shell in the latter model, which is not 
possible in the first. 

The characteristic temperature of the lattice, T*, was 
phenomenologically defined as the temperature at which 
the ratio of effective local moment n1ff{T) — Txioc{T) 



and the free moment 



= lim iiiffiT) equal 0.4 



lilff{T*)/^i^j^^^ = 0.4. xioc{T) is the local susceptibility 
of the ground state spin. Therefore, T* is a measure of 
the Kondo screening of the local magnetic or quadrupolar 
spin. 

We investigated the pair-susceptibility Xp in the in- 
termediate valence, the Kondo and the stable moment 
regime for band- fillings between quarter and 3/4-filling 
and COM momentum of q = and q — 7r(l, 1, 1). We al- 
ways performed a finite size scaling analysis for the largest 
eigenvalue of M, Eqn. (p9|), with respect to the matrix size: 
\m approaches a constant for 30 — 40 Matsubara frequen- 
cies. The superconductivity is driven by the low energy 
components of the two-particle vertex in contrary to the 
magnetic phase transitions |lj] where 400 — 500 Matsub- 
ara frequencies are needed in order to obtain a suscepti- 
bility invariant of the size of the matrix. A superconduct- 
ing transition was always obtained, but only for q = 0. 
Although Tc decreases with decreasing coupling constant 
g — PqJ — A/lsatrl, the normalized transition tempera- 
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ture tc = Tc/T* shown an increase due to the exponential 
sensitivity of T* on g as depicted in Fig. |^(a). T* decreases 
away from half- filling ric > 2: hence tc increases with Uc 
as shown in Fig. ||(b). The asymmetry of T* as function of 
ric simply stems from the asymmetry of energy of the two 
local doublets. This asymmetry disappears for £aa — 0. 

The odd-frequency vertex is attractive (repulsive) in 
the StCt (SsCs) sector as seen analytically in the weak 
coupling limit ^ from (|4|), ( ^ and ( ^8|) in combi- 
nation with an NCA approximation for n^^(iujn,iLUm',0) 
given by (|l|). 



5 Dynamical Mean Field Theory of the 
Superconducting phase 

We have established the phase boundary between the ho- 
mogeneou s, par amagnetic and the superconducting phase 
in section 2.2.2 . Only the one-particle and two-particle so- 
lutions of the effective site enter the calculation of the pair- 
susceptibility in paramagnetic phase. In this section, the 
modified DMFT equations for the superconducting phase 
are derived and approximately evaluated. New insight on 
the physical properties of this phase will be gained. 



5.1 Green Functions 

The Greens function is an 8 x 8 matrix in the spin and 
orbital space and formally written as 



G(«a;„,k) = 



^ \iujn,'k) -g(iw„,k) 
-g(-^c^„,k)]''S-l(^c^„,k) 



(59) 



where i7~^(icj„,k) — -~E_^^{—iuJn,ii)- E_ describes the 
electron, _ff the hole propagation in the paramagnetic phase 
if ^ = . The relation between the off-diagonal matrix el- 
ements is a consequence of the definition of the Green 
function and, therefore, is a fundamental relation inde- 
pendent of possible symmetries of the order parameter. 
^(ia;„,k) is a 4 X 4 matrix in the combined spin/channel 
space containing 16 linear independent anomalous Green 
functions. Since the Hamiltonian conserves magnetic and 
channel spin, the natural decomposition is given by 



£(icj„,k) = ig^[<P'a{iuJn,'k-)l + ^siiLOn,'i<.)g] 

X ["^0 («^n ' k)l -t- *c (iuJn , k)x] 

+^Q{iLUn, k)l*c(jt^n, k)x 

+^ {iUJn , k)i*s (i , k)^ 

+ *s(j^n,k)CT*c(«t^n,k)Tl 



(60) 



sector. The similarities to He'^ are striking: in He'^ the or- 
der parameter is a tensor which also can be written as 
a product of spin and orbital functions. Focusing on the 
individual symmetry sector of interest, the spin/channel 
matrix product ^(iwn, k)[^(— iw„, k)]^ is obtained in the 
StCt sector 

f(zL^„,k)[g(-zu;„,k)]'' = (*c(^c^„,k)*:(-^L^„,k)) (61) 

(*s(icj„,k)**(-iLj„,k))l 
= Tr [^(icJ„, k)^(~icj„, k)] 1 , 

where we have introduced the a new parameter matrices 



d{iujn , k) = *s {iijJn , k) («w„ , k) , 



d{iujn,k) = *;(jw„,k)**^(iw„,k) 



(62) 



and restricted our search to "unitary states" : ^c(s) (jt^^n, k) x 
*c(s) i~^^ri, k) = in obvious analogy to the treatment of 
He"^ [Q. Under the assumption that only the StCt-sector 
becomes superconducting, the full 8x8 Greens function 
is then given by 

a{^uJn,k) = 4 (Tr [^-i(ic^„,k)£-i(ic^„,k)] 
-4Tr \d{iu}n, k)^(-jw„, k)] ) 



/ K ^(it^„,k) g{iujn,k) 
[ [^(-^c^„,k)]''^-l(^c^„,k) 



.(63) 



Without directional coupling between spatial and spin/orbit 
degrees of freedom, the tensor dij separates as 



d{iuj„, k) = g{iujn, k)ns • n, 



T 

C J 



(64) 



where ng and ric"^ are constant unity vectors in spin and 
channel space, and an amplitude function g{iujn, k) = 
— g(— i(jj„, k). If the superconducting phase is chosen in 
such a way that g is a purely imaginary function of the 
arguments, and 5*(— iw„,k) = g(ia;„,k) holds, we obtain 
the Nambu Green function as 

-5(iLj„,k)g(jw„,k)]"^ 

g{iuJn, k) (^2=2 ["s^^cZ] ) i'i^n " -£'(«W„) - Efc)! 

We have reduced the 8x8 Nambu Green function-matrix 
to a standard 2x2 problem for the amplitude function 
g{z, k) and the diagonal self-energy E with the above as- 
sumptions. 



where r acts in channel space and a in spin space. The 5.1.1 DMFT One-Particle Self-Consistency Condition 
anomalous components are given by the product 

<?g(iw„, k)(pQ{iuJn, k) in the SsCs sector, <?g(iw„, k)#c(^i^nj k)As shown in the previous section, the 8x8 Nambu Green 
in the SsCt sector, f?g(ia;„, k)$s(«w„, k) in the StCs sec- function matrix can be reduced to a 2 x 2 matrix deter- 
tor and the tensor = (icjn, k)^^(iw„, k) in the StCt mining the dynamics of a homogeneous superconducting 
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phase. In order to derive DMFT equations with a purely 
local self-energy matrix, the anomalous self-energy has to 
be restricted to isotropic pairs, e. g. S'(k) = 1 



k 



(66) 
(67) 



where G {z) corresponds to the medium or dynamical 
mean field in which the effective impurity is embedded. 
It is related to a generalized Anderson width matrix by 



(68) 



and has normal and anomalous components describing 
quasi-particle propagation and pair-creation and annihi- 
lation, respectively. All components enter the calculation 
of the effective impurity. The self-energy S_ (z) is deter- 
mined via 



S = T 



1 + GT 



(69) 



a matrix generalization of the usual DMFT self-consistency 
equations, where the diagonal elements of the T-matrix 
are given by the local quasi-particle scattering matrices 
Te = Tn{z) = V^Gfiz) and n = ^22(2) = -V^Gf{-z). 
The anomalous contribution is calculated from 



Ts = TuiiuJn) = -y"^^ nf{iuJn, iuJm]0)f{iuJm) 



/3_ 

The negative sign takes into account that 

A^2{z) = -V^Gi2{z)^~V^j{z) 



(70) 
(71) 



'^mjA — 



+ 



4^ 



Fig. 6. Diagrammatic representation of the generating func- 
tional ^^rcA for the NCA in the superconduction (SC) phase: 
the first diagram is the normal contribution, the second one 
is the anomalous contribution. The wiggly lines represent the 
Green function Pa of the quadrupolar doublet, the dashed lines 
Pa for the magnetic doublet and the solid lines are the appro- 
priate components of the media Green function-matrix G_. The 
second term generates contributions to the self-energies of the 
ionic Green functions and the anomalous T-matrix. 



aa — iLO„ 



Fig. 7. Diagrammatic representation of the local anomalous 
T-matrix in NCA plotted on an imaginary time axis, from 
bottom to top. It was obtained by contracting the incoming 
lines of the diagram shown in Fig. ^. The double dashed line 
represents anomalous Cooper-pair bath function f{z) which 
obeys f{z) = -f{-z). 



5.2 Connection between the Pair-Susceptibility and 
Green Function in the Superconducting Phase 

The DMFT equations of the previous section are inde- 
pendent of the algorithm to obtain spectral information 
on the effective site. If g{z) — > holds for all frequencies 
when T Tc, the divergence of the pair-susceptibility, or 
finding an eigenvalue of A = 1 for the matrix M, Eqn. 
(Eq), is equivalent to solving (pq) and (p^) and 



(72) 



is obtained (|A = 1) oc g{iujn)) as shown in Appendix 
1^. As expected in a mean-field theory, the quasi-particle 
self-energies are not renormalized. 

However, all eigenvectors with A = 1 which le ad t o su- 
perconducting phase transitions reported in Sec. |4.3| show 
a XjiiOn behaviour in leading order reflecting the analytic 
properties of the irreducible vertex , Eqn. (^) . There- 
fore, the assumption of a simple second order phase tran- 
sition with 1^(2:) I does not hold, indicating that also 
the quasi-particle properties are strongly renormalized by 



enteringthe superconducting phase. The linearization of 
Eqns. (|6^69) with respect to g(z) is not possible below 
Tc. This implies that Heid's theorem which states 
that odd-frequency superconducting state is associated 
with a maximum in the free energy, does not apply to the 
triplet/triplet superconducting state. Its its derivation re- 
lies on the analyticity of ^(z) on the real axis. In general, 
H reads 



1 (T 



Ts - fiT.n - T/) n + e{T,Th - T 
A=il + iT, + fTs){l + hTh + fTs) 

-{eTs + fTh){m+~fT,), (74) 

which implies the renormalization of the quasi-particle 
self-energy S_^^ by the terms —hT'^ even if the we assume 
that the quasi-particle scattering t-matrices Th and are 
not renormalized close to T^. 
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5.3 Modified Non-Crossing-Approximation for the 
Superconducting Phase 

Ts is generated by the second term of the generating func- 
tional for the effective site in a superconducting medium 
depicted in Fig. ^. The diagrammatic representation of 
Ts is shown in Fig. 0. We have evaluated (^) in NCA 
by four integrals along the branch-cuts by convoluting the 
Cooper-pair media Green function f{z) with nf(iujn, i^m] 0) 
Eqn. (UW. It was shown that the instability in xp brings 
forth a solution of ( |67| - |70| ) when expanded linearly in 
the anomalous Green functions and medium (see also Ap- 
pendix 

Since we cannot use a linearized version of the DMFT 
close to Tc as discussed in the previous section, we iterated 
the full DMFT for the superconducting phase, Eqns - 
|70| ), without any further approximations for 1~T/Tc <C 1. 
The generating functional of the effective site ^ffcA 
shown in Fig. O. For T > Tc, the equations are identical 
to the normal state DMFT since g{z) iterates to zero. 



5.4 One-particle Spectra in the Triplet- Triplet Sector 




0.5 



I ' ' ' ' ' ' ' ' ' 1 

-0.05 0.05 

CO 

Fig. 9. Comparison of the imaginary part of the quasi-particle 
self-energy in the normal state, which is also used in ap- 
proximation (i), close to Tc (solid line) and the self-energy 
QmEiiiuj-iS) obtain through (H) and the full DMFT(NCA) 
(dashed line). Parameters as in Fig. 0. 
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Fig. 8. Quasi-particle media spectral function pe{i-^) ~ 
SJm e{u!~i5)/TT (solid line) and anomalous media spectral func- 
tion pf{Lo) = Qmf{Lj-i5)/n (dashed line) in DMFT(NCA) in 
the superconducting phase close to Tc. Saa ~ —2. 



We obtained a self-consistent solution of DMFT(NCA) 
equations for T close to Tc in the superconducting phase. 
For T < Tc, the ionic propagators |^,^|3^] aquire an ad- 
ditional self-energy contribution which originates from the 
anomalous conduction-electron bath / and is generated 
by the second diagram of the generating functional 'PfJ^'j^ 
depicted in Fig. H. These corrections are self-consistently 
taken into account in our numerics and modify the quasi- 
particle t-niatrices and generate an anomalous i-matrix 
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Fig. 10. Spectral functions of the quasi-particle and the 
anomalous Green function in DMFT(NCA) in the supercon- 
ducting phase for 1 — T /Tc <C 1 in the vicinity of oj = 0. The 
dotted curve shows a fit with a * Parameters as in Fig. M. 



as depicted in Fig. ^ We start the iteration with a finite 
anomalous media f(){z) — A/{z + iT) where the weight 
A is of the order of 10~^. The temperature at which a 
finite anomalous t-matrix Ts is stabilized coincides with 
the instability of the pair-susceptibility. Above Tc, f{z) 
and Tsiz) approach zero rather rapidly while iterating 
Eqns. (|l-|7g). 

The spectral functions of the effective quasi-particle 
and anomalous media is plotted in Fig. ^: a delta peak at 
a; = is growing on top of the regular part similar to the 
reported media in the single-channel Anderson lattice at 
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T = 0, shown in Fig. 3 of Ref. @. This delta peak should 
lead to a Fermi-liquid form of the quasi-particle self-energy 
Sii for T ^ 0. We clearly see the decrease of the quasi- 
particle scattering rate slightly below Tc as shown in Fig. 
^. This suppression of the quasi-particle scattering rate 
indicates the removal of the residual entropy of the nor- 
mal state: the magnetic scattering is strongly reduced due 
to forming a composite order parameter entangling local 
magnetic (quadrupolar) degrees of freedom with spin sin- 
glet /quadrupolar triplet Cooper pairs. 

The local quasi-particle and anomalous Green function 
is shown in figure |lC|. It develops a pseudo-gap for very low 
frequencies which behaves like as indicated by the dot- 
ted fit curve in Fig. |lO[ This approximate | a; | -dependency 
of the spectral functions can be understood analytically. 
It is apparent from the evaluation of Eqn. (66) that the 
frequency dependency of the spectral function is deter- 
mined by the imaginary part of the inverse of the de- 
terminant \J {z — {'^{z) -t- z))/2)^ — g^iz). Let us as- 
sume that the leading contribution to g[z) can be approx- 
imated by a temperature-broadened Lorentzian / [z + 
z sign(3mz)aT), which reflects the 1/z divergency of 77/(iw, 
at T = (see Eqn. (|2l])); the weight A is proportional to 
the order parameter, and a a dimensionless constant of 
the order of 0(1). Due to this 1/z divergency, its analytic 
behaviour has similarities to the odd frequency mean-field 
solution of the anomalous self-energy by Coleman et. al. 
pa, who found that g{z) = S12 cx I/2:. 



6 Conclusion 



finite value at the transition, the phase transition would be 
first order and the associated latent heat L = Tc{Sn — Sg) 
would be related to the change of entropy. However, in the 
presence of an anomalous medium, the quasi-particle self- 
energy is also modified and, therefore, the feedback into 
the effective site must be calculated self-consistently. The 
generating functional for the local site in a superconduct- 
ing medium is shown in Fig. |^ for the NCA. We expect 
that the self-consistent solution for the one-particle Green 
function in the superconducting phase will show a vanish- 
ing scattering rate for the quasi-particle at lim r"™ — > 0. 

In this case, a non-analytic g{z), even with an infinitesimal 
g, will be able to always produce a gap. 

It turns out, however, that the DMFT(NCA) equa- 
tions for the superconducting phase are numerically un- 
stable for T <C Tc. We believe that the instability of the 
equations is related to the known inability of the NCA to 
reproduce the Fermi-liquid phase of the periodic Ander- 
son model [0|l^ . It is clearly related to the occurrence of 
an additional very sharp delta peak in the media spectra 
(Fig.Sv Nevertheless, the solutions obtained at Tc show 
'a^'Sduction of the local free energy in the presence of a 
finite anomalous Green function and a tendency to reduce 
the scattering rate for the quasi-particles. Therefore, even 
though we cannot derive conclusively the type of the phase 
transition, we believe that there are strong hints towards 
a second order transition in the StCt sector. This phase 
transition should be associated with a minimum of the 
free energy since an energy gap develops in the spectrum 
of the local band Green function indicating an energy gain 
by condensation. 



There is a controversial debate (see and reference 
therein) as to whether an odd-frequency solution with a 
COM momentum Q = is connected to a minimum or 
a maximum of the free energy and, hence, may be ther- 
modynamically unstable. Assuming a Fermi-liquid normal 
phase and the occurrence of an infinitesimal odd-frequency 
anomalous self-energy g{z) at Tc, which implies a second 
order phase transition. Held showed that one indeed ob- 
tains an increase of the free energy in the superconducting 
phase for Q = |^ . We will give two arguments why we 
believe, in our case on the contrary, that the supercon- 
ducting phase is thermodynamically stable, even though 
we cannot prove unambiguously whether the phase tran- 
sition is first order or continuous. It is straightforward to 
extend Heid's theorem onto a non-Fermi liquid normal 
phase. An infinitesimal and analytic odd-frequency g{z) 
does lead to an enhancement of the quasi-particle DOS p 
at a; = and hence to a maximum in the free energy. How- 
ever, based on the self-consistent solution for the anoma- 
lous Green functions ( |67[]70| ), g{z) converges to a finite 
value such that T^ is approximately Fg ~ Fqp generating 
a gap in the lattice Green functions similar to the one re- 
ported earlier ||4^ . g{z) might be approximated in leading 
order by a Lorentzian with a finite weight g « ^ Fq^aT . 
The assumption of an infinitesimal value of g(z) for all 
frequencies is violated and hence Heid's theorem is not ap- 
plicable in our case. Since the order parameter then has a 



0.025 




0.5 1 1.5 2 2.5 3 



Fig. 11. The phase diagram of the two-channel periodic 
Anderson model for ric — 2.2: the superconducting transi- 
tion temperature Tc, the antiferromagnetic transition temper- 
ature Tm{AF), the spin-density wave transition temperature 
Tm{S DW) are plotted versus the effective dimensionless cou- 
pling constant g = 7rzi/|eaa|. Spin/channel triplet supercon- 
ductivity dominates the IV regime and the cross-over region 
to the Kondo regime. 
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6.1 Phase Diagram and Experimental Relevance 

The phase diagram shown for ric — 2.2 in Fig. |l^ sum- 
marizes our investigation of superconductivity and mag- 
netism in the two-channel Anderson model. Superconduc- 
tivity dominates the intermediate valence regime and the 
corresponding order parameter has StCt symmetry. An 
instability in the SsCs reported in the two-channel Kondo 
lattice could be reproduced but occurs at much lower tem- 
peratures. 

The magnetic or quadrupolar phases of the model was 
discussed in a previous publication jl^ : Fcrromagnctism 
is only found for low band fillings ric ~ 1. Spin-density 
wave phase transitions take over in the Kondo-regime at 
g w 1.3 corresponding to Saa — —2.4. The SDW wave- 
vector is continuously shifted towards nearest-neighbour 
antiferromagnetism, which is suppressed for g —^ 0. Since 
all calculations were performed in the paramagnetic phase 
of the model, the highest transition temperature defines 
the nature of the incipient order. It cannot be ruled out, 
however, that an SDW phase is replaced by, or even coex- 
ists with, a superconducting phase as found in some Ura- 
nium based Heavy Fermion compounds. Since the quasi- 
particle spectral function vanished linearly for uj 0, 
p{w) (X as depicted in Fig. |l^, for a spatially isotropic 
gap function, the specific heat obeys a power law T" below 
Tc- Additional anisotropic components oc contribute in 
the intermediate valence regime. Therefore, a exponent 
2 < a < 3 is expected: a value of a ~ 2.8 was experimen- 
tally observed on UBcia. 

In the traditional approaches to Heavy Fermion super- 
conductivity, one solely focuses on possible k-dependency 
of the gap function with nodes on the Fermi surface in 
order to explain the power-law dependency of the specific 
heat, for instance. Here, we presented a pairing mecha- 
nism which leads to a gap function with nodes in frequency 
space additionally to a k dependence (x ek. The pairing 
interaction correlates spatially extended, nearest neighbor 
Cooper-pairs to local spin or quadrupolar moments. It was 
not necessary to introduce finite center of mass momen- 
tum Cooper-pairs, i. e. Q 7^ 0, to observe the phase tran- 
sition into this state. Finite Q Cooper-pairs would imply 
a persisting current, and phases of this kind have sofar not 
been observed experimentally. 

Assuming a quadrupolar ground state of the Uranium 
ion. Fig. |ll| will show the phase diagram for quadrupolar 
order since spin and channel indices are just exchanged. 
The superconducting order parameter, however, i s no t af- 
fected by this transformation as discussed in Sec. 3.3. The 
tensor order parameter Oij, (|4^), transforms as 5*0(3) (g) 
5*0(3). In a cubic crystalline field environment, this order 
parameter reduces to Affi-Hs© Affi-^s and the phenomeno- 
logical approach by Sigrist and Rice |Q is still applicable 
in our case. 

Complicated phase diagrams as found in Thorium doped 
Ui-xThxBeis are explained by order parameters of differ- 
ent irreducible representation. Recently, the second phase 
transition of Ui_a;Tha;Bei3 was interpreted as a magnetic 
phase Q coexisting with the superconducting phase. If 
our model has any relevance for Ui-ajTh^jBeia, we would 




Fig. 12. The transition temperature Tc as function of the level 
difference Sacr and the conduction band filling. Uc = 2 corre- 
sponds to half-filling. Note the inversion of the axis scale in or- 
der to visualize more clearly the dependency of Tc for Sacr ~ 1 
and for large band fillings ric, a regime which might be of in- 
terest with respect to UBeia. 



expect an antiquadrupolar order or a quadrupolar density 
wave, since a magnetic ground state of Uranium contra- 
dicts the weak magnetic response of van-Vleck type . 

The superconducting transition temperature Tc{eaa-, n-c) 
which is shown normalized to T* in Fig. ^ is only sym- 
metric with respect to £0,0- at half filling. We have sum- 
marized all the information on Tc{s aa , nc) / A in a 3d-plot 
displayed in Fig. ^ We note that for £0,0- « 1, i. e. an in- 
termediate valence regime with quadrupolar ground state 
on the local shell as in Uranium in cubic environment, Tc 
decreases with increasing hybridization and band filling. 
A similar behaviour is observed in Vi-xThxJieis where Tc 
decreases under applied pressure and also under Thorium 
doping which is interpreted as increasing the number of 
band-electrons Q. 

The skutterudite materials PrFe4Pi2 and PrOs4Sbi2 
are also prominent candidates for the SU{2) ® SU{2) lat- 
tice model, since Pr^+ is hosted in a cubic crystalline en- 
vironment. The anti-quadrupolar ordered ground state is 
observed in PrFe4Pi2, while PrOs4Sbi2 shows supercon- 
ductivity below 2K. We expect that PrFe4Pi2 becomes 
superconducting under pressure. 
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A Composite Order Parameter 

The composite order parameter can be related to the time 
derivative of the isotropic and extended s-wave anomalous 
Green functions in the appropriate symmetry sector - see 
table |l| on page ^. The analytically continued Matsubara 
Green function for fermionic operators A and B obeys the 
equations of motion 



z GaA^) ^ {{A.B}) - G[ha],b{z) 



(75) 
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The commutator of a conduction electron with the local 
quantum numbers aa at site v, 

[H, Ci,aa] = -'ip^ao- - VaX'l^ ,^ , (76) 

is derived from the Hamiltonian (|l|), where 
c - — V e' 



ikR„ 



J / , ^ c^kacrL-ko 



(77) 
(78) 



has been used. Using the equation of motion 



(79) 

to calculate the equal time expectation value of the time 
derivative {^Cyoi<y{T)cya'<j') in the form 



= ^ E e'^""'*^nG.,_, [zujn) (80) 



yields 



The first term on the r.h.s. represents the averaged ki- 
netic energy, which is zero in the case of an odd-frequency 
anomalous Green function Gc-^^^,^,,c^i^aa{z)i while the sec- 
ond describes the coupling to the localized /-shell. Using 
the commutator 

a' 

^ ^ Vg^Q, — g^Cjyc^^fj , (82) 



and 



(Gxs,,,c„,„,„,(*^„)+ (83) 



we derive 

2VBZ:i,y) = C::: -2{E,- E^){XI,c,.„,,,) {M) 

via equation of motion and, therefore, with AE = E„ — Ea 
and Va = aV 



\AE\D'^^, -sigiv{AE)a 



= ^aVhign(AE)B^:, {i 



(85) 



The quantity B^^," (v) is defined as 

a" 

+ '^a" {X_a-a"Cua"<yCv'a'a') ■ (86) 



A relation for C"", can be derived using ([75|) and ( |82| ) 



(87) 

and the assumption that the anomalous Green function 
is odd in frequency. If only pairs with COM momentum 
Q = condensate, the condition 

holds. Then f^^",' is equal to T^^'^', defined in Eqn. (|l]). 
It is straightforward to show that 

((S P*'^)-(r P^^*))=-^ai?r. 



and to obtain the equation ( [42| ) from (g5|) and (|88|): 
l/i-Bl ^ sign(Z\i;) 



y2 



■5]aaT,%"(89) 



Eqns. ( p5| ) and (|8^) in conjunction also yield Eqn. ( p4| ) on 
page H for the order parameters in section IV. 



B Calculation of the Odd- Frequency 
Pair-susceptibility 

Assuming we know the irreducible particle-particle ver- 
tex iZ^^^ for the appropriate symmetry and only consider 
isotropic pairs, then the pair-susceptibility of the compos- 
ite order parameter is proportional to 



XP 



x(q,o) 



1 



i-^£ff JO)x(q,0) 



(90) 

Since x(q, ii^n) is diagonal in Matsubara frequency space. 



x(q, 0) is a well defined quantity so that a new symmet- 
rical matrix M 



(91) 



can be introduced. It is diagonalized to obtain the base 
I A), Eqn (H), M|A) = A|A) and to rewrite xp(q) 

Xp(q) = i(a;„|y^^(q,0) [i - M] y^^(q, 0)|u;„) , 

(92) 

where |w„) = (• • • , a;„, • • ■)'^ . This eigenvalue equation is 
viewed as the analog to the linearized Eliashberg equation 
in the standard theory of superconductivity, even though 
no reference to Migdal's theorem was made here. We can 
use all eigenvectors |A) as a complete base in frequency 
space 

i = ^|A)(A| (93) 

A 
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and insert it into 193) to obtain 



(94) 



Since ^x(Qi i^7i)\^n) is an odd-frequency vector, only odd- 
frequency eigenvectors |A) contribute to Xp(q)- The lead- 
ing contribution is given by the eigenvector |Ac) for which 
|1 — A| has a minimum. Hence, we can write 



xp(q) = icl 




regular terms 



(95) 



It is now obvious that the pair-susceptibihty Xp(q) di- 
verges when Ac ^ 1, and that Ps(q) is positive in the 
paramagnetic phase. 



C One-Particle Self-Consistency Condition for 
Local Cooper Pairs 

In DMFT, a purely local self-energy leads to a self-consi- 
stency condition (SCC), which states that the k-summed 
one-particle lattice Green function equals the Green func- 
tion calculated for an effective site with excitations into a 
surrounding effective medium: 



k L 

;G(z)]"'=[g(z)]"+^^(z). 



1 -1 



(96) 
(97) 



G°(k, z) denotes the Green function in the absence of any 
interaction, and G^{z) denotes the medium or dynamical 
mean field which enters the effective site problem. It is 
related to the Anderson width matrix by 



(98) 



where is the 8x8 generalization of (72 appropriate 
to the description of the superconducting state with spin 
and channel degrees of freedom. It was assumed in Sec. |3.1| 
that the orientation in spin and channel space is described 
by the two unitary vectors. Hence, we use a 2 x 2 matrix 
formalism for the dynamics. Close to Tc, e.g. for g{z) — > 0, 
Eqn. (pa) is approximated by 



G (z,k) 



1 



{z~ ek)(^ - Sh+ Sk) 

z- Sh + e-k 9{z) 
g{z) z- Se- £-k 

S{z) ; = -Si-z) 



(100) 



in linear order in g{z). The subscripts e and h refer to elec- 
trons or holes. We immediately recognize, in conjunction 



with Eqn. (|9q), that the quasi-particle properties are not 
renormalized close to Tc- This is typical for a mean- field 
theory as the DMFT. With the local electron and hole 
Green functions 



k 

M^)-^E 



1 



TV ^ z - - ek ' 

k 



N ^ z-Sh + Ei, 

k 



(101) 

= -e(-z) , (102) 



we obtain the off-diagonal matrix element 

e{z) + h{z) 



/W- Giz) 



12 2z ~ Se — Sh' 



(103) 



which we insert into the SCC condition ( p7| ) with the re- 
sult 

/>) _ fiz) 



i{z)h{z) eiz)h{z) 



9{z) 



(104) 



where e = 
and h ~ 



Giz) 



corresponds to the electron medium. 



G (z) corresponds to the hole medium. / — 

. — J 22 

G (z) denotes the anomalous medium in the super- 

.=c J 12 

conducting phase. g{ioJn) is generated by the particle- 
particle irreducible vertex r^P(iujn,ii^m) 



g(lUJ„) ^ J '^''[lUJn,lUJm)J[it^m) (105) 

which is local in the DMFT. The anomalous Greens func- 



tion / can be replaced through the SCC (104) yielding the 
compact vector equation 



where 
and 



9 = -^L'%Mr'i' 



e{iuJn)h{iojn) 



Here, we introduced the short hand notation 
We use the irreducible vertex (p8|) and note that x 



-A to obtain 



(106) 

(107) 
(108) 

(109) 

'■loc 

(110) 



(99) It has been discussed in section 2.2.3 that the local two 



particle Green function ILi^^ ( |22| ) can be written as 



n. 



-PA+BT B, 



(111) 



where T is the local two-particle T-matrix cumulant 



(112) 
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and 11^ the local two-particle f-Green function. Hereby, 



[£]n,m = F{iujn)F{-iujn)5n,m- Hencc, Eqn. ( |105D reads 



9^ -^A'^ RZ^l- 
With the matrix equation for the vector 

[j—t - 



(113) 



(114) 



which can be derived via equation of motion or functional 
derivative, and 



1 



e 1-1- eTe 
we obtain 

9 



and 



h 
h 



1 + hTh 



(115) 



1 



(t, - TeThf) (116) 



for each frequency iuJn which is just the off-diagonal ma- 
trix element of the 2x2 self-energy matrix S_ 



S = T 



l + GT 



(117) 



in the limit f,f^O where the one-particle T-matrix is 
given by 

In Eqns. ( |115| - |118| ), the energy argument icj„ was dropped 
for simplicity. In general, S_ reads 

£ + hiT^n - T]) - ~f{T,Tn - Tf) \ 
= A\Ts- fiTcTh - T^) n + i{T,n - T^) } ' 
A={l + eT, + /T,)(l + hTu + ~fT,) 

-{eTs + fn){hT, + fT,) . (120) 

As mentioned above, g{(\, i^n) is generated by the par- 
ticle-particle irreducible vertex rPP{iiVn, «w„/), Eqn. (105) 



g(q,iLj„) = --^ ^ rPP{iu;n,ii^7i')h^ci{ii^n') (121) 

which is local in the DMFT. On the other hand, we can ex- 
pand the k-dependent anomalous Green's function /k,q(it^n) 
close to Tc up to linear order 

/k,q(«w„) ^ -Gi^+q{iujn)G-k{-i^n) 9{^,i^7i) + 0{g^) 

(122) 

in the anomalous self-energy (7(q, ia;„) for Coo per pairs 
with a center of mass momentum q. Inserting (|l22|) into 



121) and identifying 

Xo^(iLj„;w„,q) = — '^G\^+c^{^|^n + ivn)G-\^{-iu!n) 



(123) 



yields 

g(q,zw„) = i^rPP(itj„,iw™)xo^(q,0)g(q,zcj,„) . 

(124) 

If we set g oc |eliash) = ^x(q = O^O)) ^ 1^ — 1)j where 
I A) was a solution of the eigenvalue equation (^ Ml-^) = 
A I A), the vector jeliash) is a solution of the linearized 
Eliashberg equation (fz^). Since |A = 1) corresponds to a 
diverging pair-susceptibility xp (H^); ^ singularity in xp 
is equivalent to a solution of the SCC ( |97| ) close to Tc with 
non- vanishing off-diagonal matrix elements. 
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